Pulsed Nuclear Pumping and Spin Diffusion in a Single Charged Quantum Dot 
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We report the observation of a feedback process between the nuclear spins in a single charged 
quantum dot and its trion transition, driven by a periodic sequence of optical pulses. The pulse 
sequence intersperses off-resonant ultrafast pulses for coherent electron-spin rotation and resonant 
narrow-band optical pumping. The feedback manifests as a hysteretic triangle-like pattern in the 
free-induction-decay of the single spin. We present a simple, quasi-analytic numerical model to 
describe this observation, indicating that the feedback process results from the countering effects of 
optical nuclear pumping and nuclear spin-diffusion inside the quantum dot. This effect allows dy- 
namic tuning of the electron Larmor frequency to a value determined by the pulse timing, potentially 
allowing more complex coherent control operations. 



Optically controlled quantum dots (qds) are in many 
ways similar to atomic systems, and are therefore of- 
ten regarded as strong candidates for solid-state quan- 
tum information processing. However, one key feature 
distinguishing qds in group ni-v semiconductors from 
atomic systems is the presence of a large nuclear-spin en- 
semble 1 j. Nuclear spins cause adverse effects such as 
inhomogeneous broadening and non-Markovian decoher- 
ence processes. However, nuclear spins may play useful 
roles as well. Although methods to use QD nuclear spins 
directly as a quantum memory remain challenging due 
to the difficulty of achieving sufficiently high levels of nu- 
clear polarization, nuclear spins may provide novel meth- 
ods for the dynamic tuning and locking of electron spin 
resonances for electrons trapped in qds. 

Several examples of manipulating nuclei to improve 
electron spin coherence have recently been observed. In 
electrically controlled double qds, transition processes 
between electron singlet and triplet states allow the ma- 
nipulation of interdot nuclear spin polarizations, improv- 
ing coherent control [21 El H]. In single qds under mi- 
crowave control, nuclear effects dynamically tune the 
electron spin resonance to the applied microwave fre- 
quency [5j. Tuning effects are also observed in two-color 
continuous-wave (CW) laser experiments, in which the 
appearance of coherent electronic effects such as popu- 
lation trapping are modified by nonlinear feedback pro- 
cesses with nuclear spins jH [7]. Finally, nuclear spins 
have been shown to dynamically bring ensembles of in- 
homogeneous QDs into spin-resonance with a train of ul- 
trafast pulses [3 E] ■ 

Here, we describe a related but different manifestation 
of the non-Markovian dynamics occurring between a sin- 
gle electron in a QD and the nuclear bath with which 
it interacts, with new possibilities for use in controlling 
nuclear effects. The effect occurs when measuring the 
familiar "free- induction decay" (fid) of a single spin in 
a single QD under pulsed control. The Larmor frequency 



of the electron spin is dynamically altered by the hyper- 
fine interaction with QD nuclei; the nuclear polarization 
is in turn altered by the measurement results of the FID 
experiment. The result is a feedback loop in which the 
nuclear hyperfine field stabilizes to a value determined 
by the timing of the pulse sequence. In what follows, 
we show the experimental manifestation of this feedback 
loop and present a numerical model for the effect. 

The fid or Ramsey interferometer experiment pro- 
ceeds by tipping a single electron spin from the pole of 
its Bloch sphere onto the equator with a 7r/2 rotation, 
allowing it to precess for a time r, and then tipping it 
with a second 7r/2 rotation into a state with polarization 
depending sinusoidally on r. We accomplish this exper- 
iment in a single, charged, self-assembled InAs QD using 
periodically applied pairs of ultrafast optical pulses, in- 
terspersed with optical pumping. The molecular-beam- 
epitaxially grown qds are centered in a planar micro- 
cavity with 24 and 5 pairs of GaAs/AlAs A/4 mirror 
pairs below and above the cavity, respectively. The cavity 
quality factor was roughly Q ss 200, and the single-sided 
design of the cavity directed most emission towards the 
collection optics. A layer of Si dopants placed 10 nm be- 
low the QDs served to probablistically dope the qds with 
electrons. Several prior experiments have shown that sin- 
gle ultrafast pulses achieve coherent, single-spin rotations 
in charged qds [TUl [HI H2]- The rotation occurs on pi- 
cosecond timescales, which is effectively instantaneous in 
comparison to electron-spin decoherence times or nuclear 
spin dynamics. The optical pumping in our experiment, 
previously described in Ref. HTJ is accomplished via a 
trion transition. It serves to repeatedly re-initialize the 
spin into one particular polarization state prior to the 
ultrafast pulses. Spin measurement is accomplished by 
single-photon-counting during the optical pumping step; 
the photon count-rate is proportional to the probability 
of exciting a trion state. This count-rate averaged over 
many sequential experiments is expected to vary sinu- 
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FIG. 1: (Color online) (a) Experimental Ramsey fringe count- 
rate as a function of two-pulse time delay r, measured using 
the techniques of Ref. 1111 This data was taken at a magnetic 
field of 4 T, but similar effects are seen at different magnetic 
fields, (b) Average electron polarization as a result of the peri- 
odic pulse sequence used to generate this data. Optical pump- 
ing increases the polarization for a duration T — 26 ns. The 
saturation polarization, which would be reached for T — > oo 
(dashed line) is S^; in time T only polarization Sf is reached. 
After pumping and a short delay, a picosecond pulse indicated 
by a green (grey) arrow nearly instantaneously rotates the 
electron spin to the equator of the Bloch sphere ({S z ) = 0); a 
time t later a second pulse rotates the spin to achieve electron 
polarization Sf , depending on the amount of Larmor preces- 
sion between the pulses. The theoretical count-rate C(u,r) 
of Eq. |l| is found as Sf — Sf in steady-state conditions, (c) 
Experimental Ramsey fringe count-rate as r is continuously 
scanned forward and then backward, showing clear hysteresis. 

soidally with r, resulting in "Ramsey fringes," and to 
decay on a timescale called T% due to dephasing or deco- 
herence effects. 

In contrast to earlier work [TT] , the resonant laser used 
to optically initialize the single spin is gated off during 
the free precession of the electron spin using an electro- 
optic modulator. This eliminates the previously observed 
source of decay of Ramsey fringes, leading to an expected 
longer T 2 * decay. In particular, a Markovian picture of 
nuclear spin noise processes predicts a Gaussian decay 
with T 2 * timescale on the order of a few nanoseconds. 

However, such a Gaussian decay was not observed. 
Figure [T| shows the result of the fid experiment. The 
top three traces show the fringes seen as r is scanned in 
the forward direction (increasing the delay r by moving 
a retroreflector outwards with a scanning stage) , and the 
bottom three correspond to scanning in the backward di- 
rection (reducing the delay r by moving the retroreflector 
inwards). The oscillatory fringes, rather than decaying, 
evolve into a sawtooth pattern at high values of r [15] . 
and show hysteresis depending on the direction in which 



t is changed. The break between traces in Fig. [l|a) is 
due to the need to manually move the scanning stage. 

This data is the result of two competing processes: 
changes in the average nuclear hyperfine shift lo due to 
trion emission, in conjunction with the motion of that 
magnetization due to spin diffusion. In what follows, we 
first qualitatively describe these physical processes and 
explain how they lead to our data, and then we present 
equations to formally model the dynamics quantitatively. 

One important assumption is a separation of dynamics 
into three very distinct timcscalcs. The fastest timescale 
is the pulse sequence and resulting electron-spin dynam- 
ics, repeated continuously with a repetition period of 
143 ns, shown in Fig. [ljb). This is much faster than 
the nuclear dynamics we consider, which are presumed 
to occur on millisecond timescales. Finally, the averag- 
ing timescale of the measurement is much longer still, on 
the order of several seconds, allowing the nuclei substan- 
tial time to reach quasi-equilibrium. 

Processes that change the total nuclear magnetization 
at the high magnetic fields used here (4 to 10 T) are 
unlikely to be due to the flip-flop terms of the contact 
hyperfine interaction of the ground-state electron in the 
QD, as its energy levels are known to be narrow (on the 
order of Ti/T 2 , with T 2 ~ 3 /xs) leaving few viable path- 
ways for energy-conserving nuclear-spin flips. In con- 
trast, the dipolar interaction between a trion's unpaired 
hole and a nuclear spin may induce a spin-flip with the 
nuclear Zeeman energy compensated by the broad width 
of the emitted photon (j/2tv ~ 0.1 GHz). Fermi's golden 
rule allows an estimate of the rate at which a trion hole 
(at position rh, with gyromagnetic ratio <?h) polarized 
along the sample growth axis (orthogonal to the mag- 
netic field) randomly flips a nuclear spin at position r 
in a spatially flat QD during spontaneous emission, with 
the photon energy density of states negligibly changed by 
the Zeeman energy of the nucleus. The result is T(r) sa 
(9 At 2/1287r)( AtB5h /B ) 2 7 (|r-r h |- 3 ) 2 ~ 1/(20 ms), where 
the brackets refer to an average over the hole wavefunc- 
tion. Nuclear polarization due to this process has been 
considered before in the modeling of similar effects 6 , 9 . 
The change in the nuclear magnetization induced by the 
trion is random; spins may flip in either direction, leading 
to a random walk. However, the rate of this walk is it- 
self modified by the probability that the trion is excited. 
This probability is a function of the pulse sequence; for 
the sequence shown in Fig.[TJb), we take it as 

[1 - exp(-/3MT)] [1 - cos[^o + lj)t}} 
1 ' 1 p 1 - cos[(w + u) T ] exp(-/?(w)T) ' 

(1) 

Here, loq is the Larmor frequency of the electron spin in 
the absence of nuclear shifts, lo is the randomly drifting 
Overhauser shift, /3(w) is the rate of optical pumping, 
T is the pumping time, and S p is the saturation value 
of the polarized spin, equal to 1/2 for perfect pumping. 
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FIG. 2: (Color online) Count-rate C(lu,t) as a function of 
Overhauser shift u and two-pulse delay r. The lighter grey or 
greener areas indicate where a higher count-rate is expected. 
Oscillations in the horizontal directions at frequency u>o + u 
are due to Ramsey interference; the Gaussian envelope in the 
vertical direction is due to the reduction of optical pumping 
with detuning. Superimposed in black is a line indicating 
where dui/dt = according to Eq. Q. Superimposed on this 
line are the solutions to this equation which result as r is 
scanned forward (yellow or light grey) and backward (white). 



This function is plotted in Fig. [2] C(lj, t) oscillates sinu- 
soidally with increasing r due to the spin's Larmor pre- 
cession; the Overhauser shift lj affects the frequency of 
Larmor precession. For large values of |w|, the trion tran- 
sition shifts away from resonance with the optical pump- 
ing laser, leading to reduced pumping efficiency and trion 
creation. If magnetization drift due to trions were the 
only process, u> would randomly drift to whatever value 
is needed to null C(lu,t), thereby slowing and stopping 
the drift. As a result, no fringes would be observed. 

The second process which counters this drift is the 
presence of nuclear spin diffusion. When trion emission 
pushes u> to too large a value, nuclear dipolar interactions 
"flatten" the nuclear magnetization. As a result, the shift 
to is "pulled" back to a low value, countering the tendency 
of trion emission to push oj away from zero. The stable 
quasi-equilibrium value of lo resulting from the balance 
of these processes lives on the edge of the fringes shown 
in Fig. [2] the nuclear polarization "surfs" along the edge 
of this function as r is changed. As r is increased, |a>| in- 
creases causing the observable photon count to decrease 
due to the reduced degree of optical pumping. When \u>\ 
is so high that pumping is ineffective {(3{u>) — > 0) and 
the trion-induced walk stops, spin-diffusion causes the 
system to drift back to a new stable magnetization at a 
lower value of \uj\, and the process continues. 

These processes may be formally modeled by a diffu- 
sion equation for the nuclear distribution. In this model, 
the nuclear magnetization at each nuclear site j is a ran- 
dom variable, Mj. A probability distribution function 
(pdf) /(mi,TO2, ■ • • ;t) — /(m; t) gives the joint proba- 
bility that the nuclear magnetization at each position is 



Mj = rrij at time t. 



The Overhauser shift is then also a 



A{jCj) is the electron hyperfine field at the position Tj of 
nucleus j. The average value of fl at time t is written 
(Q) = u)(t) and is found by averaging over all possible 
values of each Mj, weighted by the joint PDF /(m;t). 
This joint PDF obeys the equation 
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dm 2 - 



(2) 



The first term, in which the sum over k is the sum over 
neighbors of j, describes the dissipative component of nu- 
clear spin diffusion with diffusion rates Djk- The second 
term describes the random walk of the magnetization at 
each location r due to stochastic nuclear spin-flips from 
the trion hole-spin; the constant Fj models the fluctuat- 
ing component of nuclear spin diffusion, a term needed 
to understand single QD T 2 * effects in the absence of the 
nonlinearities we consider here. 

A detail omitted from this model is the degree to 
which the nuclear spin diffusion is slowed or directed 
by the presence of the electronic hyperfine field. The 
magnetic field gradient introduced by the electron does 
not completely freeze diffusion in this case, because the 
electronic wave-function is large in comparison to the 
nuclear- nuclear spacing. Reference EH indicates methods 
by which this could be modeled. However, our work solv- 
ing Eq. ([2]) including such methods added greater com- 
plication and little insight into the observed data. 

Remarkably, the data of Fig. [T] can be understood by 
examining just the average shift uj(t) using two simpli- 
fying assumptions. In the first assumption, the electron 
wavefunction is modeled as a surface with a hard bound- 
ary, and the first term of Eq. ^ acts to decrease u>(t) by 
modeling diffusion through that boundary. This results 
in the equation 



dui 
~dt 



on- 



■[SlC(Cl,t)] 



(3) 



The constant k depends on the electronic wavefunction 
and the rate of nuclear diffusion, but we treat this pa- 
rameter as adjustable rather than attempting a micro- 
scopic description. The constant a is formally given by 

Unfortunately, Eq. ^ is not a closed system of equa- 
tions, because it still requires full knowledge of /(m; t) 
to solve. However, if C(f2, r) is a sufficiently flat function 
of f2 in comparison to the width of /(m; t), then we may 
treat C(£1,t) as roughly constant at C(u>(£),t) over the 
small width of f(m;t). This is our second assumption, 
leaving the final equation 



random variable f2, defined by f2 = A{rj)Mj, where 



d 2 

-kuj + a— —\uiC(lj, f) . 

OUJ z 



(4) 
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FIG. 3: (Color online) The modeled (a) countrate or Ramsey 
amplitude C(uif, r), (b) Overhauser shift cjf, and (c) Pumping 
strength f3{ui). The dotted line in (a) is the expected Ramsey 
fringe in the absence of nuclear effects. The traces in (b) are 
the same as those in Fig. [2] The blue or darker grey line 
corresponds to scanning r forward, and the red or lighter 
grey line corresponds to scanning r backward. 



Invoking our separation of timcscalcs, wc presume ui{t) 
evolves from its initial value (set by the last chosen value 
of r) to a quasi-equilibrium final value u>{. This final 
value determines the expected count rate C(u){, r) at this 
value of t. We solve by assuming w(0) = at the first at- 
tempted value of r, and then we scan r up and then down 
as in the experiment, finding the steady-state solution of 
Eq. Q at each value. 

Figure [3] shows the modeled C(ujf, r), Wf, and /3(wf) as 
a function of r. This particular model used n/a = 10 4 , 
which reproduces the qualitative shape of the data quite 
well, and quantitatively reproduces the location where 
sinusoidal fringes evolve into sawtooth-like fringes. Qual- 
itative differences are dominated by the random condi- 
tions that develop when the stage is moved on its rail, 
forming the breaks between data sets in Fig. [T|a) . Details 
of the shape of the waveform are related to the assumed 
form of the optical absorption. For simplicity, we have 
used (3{lj) = /3 exp(~w 2 /2CT 2 ), with a/2ir = 1.6 GHz 
and Po — 3/7 1 for known pumping time T = 26 ns, which 
roughly matches the experimentally observed count-rate 
when scanning the pump laser across the resonance. The 
real absorption shape is difficult to observe directly since 
hysteretic nuclear pumping effects also appear in absorp- 
tion experiments with scanning CW lasers, as reported 
elsewhere [5J[7]. 

This effect may be useful for future coherent technolo- 
gies employing QDs. This pulse sequence may serve as a 



"preparation step" for a qubit to be used in a quantum in- 
formation processor, as it tunes the qubit to a master os- 
cillator [14] and narrows the random nuclear distribution, 
assisting more complex coherent control [3J [31 IH HJ [HI IH] • 
In particular, the ability to control a single electron 
with effectively <5-function-like rotation pulses introduces 
strong potential for dynamical decoupling |151 116j . but 
many schemes, especially those that compensate for pulse 
errors such as the Carr-Purcell-Meiboom-Gill (CPMG) 
sequence, require some method to tune the Qd's Larmor 
period to an appropriate division of the pulse-separation 
time. 

In conclusion, we have observed nonlinear nuclear feed- 
back effects in a single charged QD resulting from the 
countering processes of random nuclear walks driven by 
trion creation, the finite width of optical absorption, 
Overhauser-shifted Larmor precession, and nuclear spin 
diffusion. This feedback may be employed for tuning the 
electron Larmor period to a particular pulse separation 
time for more complex pulses sequences. Although the 
model we have presented is simple and replicates the data 
well, it is insufficient to describe the processes in a QD 
under all possible pulse sequences. In particular, these 
nonlinear effects arc highly suppressed in spin-echo mea- 
surements |17j . even though trion creation follows a sim- 
ilar nonlinear function to Eq. (JT|). Future work will in- 
volve extending this model to explain the non-Markovian 
effects of nuclei under more complex pulse sequences, as 
well as exploiting it to extend QD-based quantum mem- 
ories. 
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